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ABSTRACT 



Aims. One of the puzzling findings in the study of the chemical evolution of (post-)AGB stars is why very similar stars (in terms of 
metallicity, spectral type, infrared properties, etc. . . ) show a very different photospheric composition. We aim at extending the still 
limited sample of s-process enriched post-AGB stars, in order to obtain a statistically large enough sample that allows us to formulate 
conclusions concerning the 3rd dredge-up occurrence. 

Methods. We selected two post-AGB stars on the basis of IR colours indicative of a past history of heavy mass loss: IRAS 08281- 
4850 and IRAS 14325-6428. They are cool sources in the locus of the Planetary Nebulae (PNe) in the IRAS colour-colour diagram. 
Abundances of both objects were derived for the first time on the basis of high-quality UVES and EMMI spectra, using a critically 
compiled line list with accurate log(gf) values, together with the latest Kurucz model atmospheres. 

Results. Both objects have very similar spectroscopically defined effective temperatures of 7750 - 8000 K. They are strongly carbon 
and s-process enriched, with a C/O ratio of 1.9 and 1.6, and an [ls/Fe] of +1.7 and +1.2, for IRAS 08281-4850 and IRAS 14325-6428 
respectively. Moreover, the spectral energy distributions (SEDs) point to heavy mass-loss during the preceding AGB phase. 
Conclusions. IRAS 08281-4850 and IRAS 14325-6428 are prototypical post-AGB objects in the sense that they show strong post 3rd 
dredge-up chemical enrichments. The neutron irradiation has been extremely efficient, despite the only mild sub-solar metallicity. This 
is not conform with the recent chemical models. The existence of very similar post-AGB stars without any enrichment emphasizes 
our poor knowledge of the details of the AGB nucleosynthesis and dredge-up phenomena. We call for a very systematic chemical 
study of all cool sources in the PN region of the IRAS colour-colour diagram. 
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1. Introduction 

Post-AGB stars are key objects in the study of the dramatic 
chemical and morphological changes of objects on their ascent 
on the Asymptotic Giant Branch (AGB) and subsequent evolu- 
tion. In this paper we report on our ongoing research to study 
the AGB chemical evolution by a systematic study of post- 
AGB photospheres. Spectra of post-AGB stars are much easier 
to study than their AGB precursors for several reasons. First, 
their atmospheres do not show the large amplitude pulsations 
and the large mass loss rates that characterise AGB atmospheres. 
Second, their photospheres are hotter, so atomic transitions pre- 
vail in post-AGB spectra, while molecular transitions prevail in 
AGB star spectra. This allows us to quantify the chemical con- 
tent of a very wide range of trace elements. Unfortunately, post- 
AGB stars evolve on a very fast track and whether the current 
Galactic sample is representative is one of the most difficult is- 
sues when interpreting the results of single post-AGB stars in 
the broader context of stellar evolution. Therefore, extending the 
still limited sample of post-AGB stars is indispensable in order 
to attain a statistically large enough sample that allows us to for- 
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mulate conclusions concerning the (post-)AGB evolution in gen- 
eral. 

During the past decennia, it has been realised that post-AGB 
stars are chemically much more diverse than previously thought. 
Binary objects tend to have a totally different photospheric com- 
position than single objects, showing some degree of depletion 
of refractory elements in their photosphere (see e.g. Maas et al. 
2005). The single objects, on their side, are also far from chem- 
ically homogeneous. Some objects are the most s-process en- 
riched objects known to date (e.g. Reyniers et al. 2004, and ref- 
erences therein), while others are not enriched at all. This di- 
chotomy is very strict, in the sense that mildly enhanced objects 
do not exist, except for a few rather atypical objects. Chemical 
evolutionary AGB models do predict that there is a minimum ini- 
tial mass for the 3rd dredge-up to occur, at around 1.4M (e.g. 
Straniero et al. 2003). Hence post-AGB stars without any helium 
burning products in their photosphere are theoretically expected. 
However a strict dichotomy is not expected, since a more grad- 
ual transition between non-enriched and enriched objects is pre- 
dicted if the transition from an O-rich AGB to a C-rich AGB star 
occurs over many thermal pulses. Furthermore, the s-process en- 
riched sub-class of post-AGB objects imposes another unsolved 
problem. These stars exhibit a large spread in s-process effi- 
ciency, as they do not obey the expected anti-correlation between 
metallicity and s-process efficiency (see Fig. 4 of Van Winckel 
2003). In other words, examples exist of very similar post-AGB 
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Table 1. Basic parameters of the two objects discussed in this study. 



IRAS 08281-4850 14325-6428 



Equatorial coord. 


tt2000 


08 29 40.552 


14 36 34.375 




£2000 


-49 00 04.33 


-64 41 31.09 


Galactic coord. 


1 


266.08 


313.87 




b 


-5.82 


-4.08 


Vis. magnitude 


V 


14.1 


11.9 


Spectral Type 




A9I 


A8I 



Source SIMBAD, except magnitude and spectral type 



(this work). 



Table 2. Observational log for the images obtained at the Euler tele- 
scope at ESO, La Silla. 



IRAS 


Date 


filter 


exptime 


airmass 


08281-4850 


2005-10-07 


UG 


300.2 


1.473 






BG 


20.9 


1.423 






VG 


150.5 


1.408 






RG 


60.6 


1.392 






IC 


98.9 


1.351 


14325-6428 


2005-08-04 


UG 


399.3 


1.322 






BG 


299.5 


1.311 






VG 


241.6 


1.292 



stars (in metallicity, spectral type, infrared excess, etc. . . ), but 
with a totally different photospheric abundance pattern (see for 
example Fig. 1 in Van Winckel & Reyniers 2001). This result 
dramatically illustrates that the 3rd dredge-up phenomenon is 
not yet fully understood. 

In this paper, we selected two objects with far infrared 
colours typical of PNe from the IRAS point source catalogue. 
Apart from PNe, only post-AGB stars are typically found in 
this part of the colour-colour diagram (van Hoof et al. 1997). 
Furthermore, the selected objects were not detected in the radio 
continuum above a detection limit of 3 mJy (Van de Steene & 
Pottasch 1993, 1995). Hence they do not seem to have evolved 
to the PN stage as yet and are very likely post-AGB stars. We 
obtained JHKL images of the candidates with CASPIR on the 
2.3-m telescope at Siding Spring Observatory in Australia in or- 
der to assure the correct identification of the IRAS counterparts 
and obtain accurate positions (Van de Steene et al. 2000). Two of 
these infrared selected post-AGB stars which are presented here 
(see Table 1) have counterparts in the USNO catalogue so that 
their visual magnitudes are also known. They were sufficiently 
bright to be observed in the optical at high resolution. 

The paper is organised as follows: in the next section, we 
discuss the spectral energy distribution of the two objects, and 
quantify the total reddening towards the two sources. In Sect. 
3 we briefly discuss the observation and reduction of the high- 
resolution spectra, while in Sect. 4 we deal with the abundance 
analysis. Sect. 5 is a section devoted to the diffuse interstellar 
bands in the spectrum of IRAS 14325-6428. In the discussion 
(Sect. 6), we mainly focus on the s-process abundances. We end 
the paper with conclusions. 

2. Spectral energy distribution 

2.1. Geneva Photometry 

Images were obtained for both objects with the Swiss Euler tele- 
scope at La Silla (ESO) in 2005. The observing log is presented 
in Table 2. The objects were unresolved in all images. The im- 
ages were bias subtracted and flatfielded in iraf. The photome- 
try package was used to do the aperture photometry. In case of 
crowded fields daophot was used to subtract the neighbours. One 
to three standard stars per band per star were available for the 
photometric calibration. The resulting magnitudes are presented 
in Table 3. 

2.2. Additional Photometry 

We collected additional photometry from various surveys: the 
Tycho-2 Catalogue (H0g et al. 2000), the 3 rd data release of 
DENIS (Epchtein et al. 1994), 2MASS (Skrutskie et al. 1997), 
MSX (Egan et al. 2003), and the IRAS point source catalogue 



(Beichmann 1985). We also obtained literature data from Garcia- 
Lario et al. (1997, hereafter GL97) and Van de Steene et al. 
(2000, hereafter VdSOO). The data are presented in Table 3. 
Near-IR data of IRAS 08281-4850 were present in the 2 nd re- 
lease of the DENIS data set, but was deleted in the 3 rd release. 
The data from the 2 nd release have not been retained in the anal- 
ysis. In order to construct the spectral energy distribution (SED), 
we first converted the Tycho-2 magnitudes to the Johnson system 
using the formulas: 

B. r = B T - 0.240(B r - V r ) 

V y = V r - 0.090(B r - V r ) 

which were derived from Perryman et al. (1997). The Gunn- 
r magnitude from the Geneva system was converted to the 
Johnson system using the formula: 

R J =r-0.43-0.15(B / -V / ) 

(Kent 1985) where Bj - V/ can be calculated from the Geneva 
photometry using Eqs. 12 and 13 from Harmanec & Bozic 
(2001). After that, the photometry was converted to AF\ using 
the absolute flux calibrations given in Table 4. For the GL97 ob- 
servations we adopt the Johnson photometric calibration, while 
for the VdSOO observations we adopt the MSSO photometric 
calibration. The resulting spectral energy distribution is shown 
in Table 5. 

2.3. Fitting the SED 

The SEDs of the post-AGB stars suffer from interstellar as well 
as circumstellar extinction. We assumed that this total extinction 
can be described by an Ry = 3.1 galactic extinction curve as de- 
fined by Fitzpatrick (1999). We furthermore assumed that the in- 
trinsic spectrum of the stars can be described by an Atlas stellar 
atmosphere model (Castelli & Kurucz 2004) with the parameters 
given in Table 8. We then used an iterative procedure to deter- 
mine the amount of total extinction Ay and the bolometric flux 
L/(4ttD 2 ) by minimizing the quadratic residuals of all observa- 
tions between 0.4 and 4 fim. The results are given in Table 6. 
Plots of the Kurucz stellar atmosphere models combined with 
the dereddened observed fluxes are shown in Figs. 1 and 2. 

3. High-resolution spectra: observation and 
reduction 

High-resolution, high signal-to-noise optical spectra of the two 
program stars were taken in the framework of our ongoing pro- 
gram to study the photospheric chemical composition of stars 
in their last stages of evolution (e.g. Reyniers et al. 2004; 
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Table 3. Broadband photometry for IRAS 08281-4850 and 
IRAS 14325-6428. IRAS 14325-6428 is contained twice in the 
DENIS catalogue. The data sets do not agree within the error margins, 
which may point to variability of the source. Both data sets have been 
retained in the analysis. 



phot, band 


08281- 


4850 


14325 


6428 




mag 


Amag 


mag 


Amag 


Euler Geneva U 


17.81 


0.15 


14.58 


0.25 


Euler Geneva B 


14.76 


0.10 


12.08 


0.10 


Euler Geneva V 


14.09 


0.05 


11.89 


0.05 


Euler Gunn-r 


13.23 


0.15 






Euler I c 


12.20 


0.05 






Tycho-2 B T 






12.676 


0.209 


Tycho-2 V r 






12.065 


0.169 


DENIS Gunn-i 






10.297 


0.04 


DENIS J 






9.202 


0.07 


DENIS K s 






8.445 


0.09 


DENIS Gunn-i 






10.431 


0.02 


DENIS J 






9.292 


0.05 


DENIS K, 






8.604 


0.05 


2MASS J 


10.603 


0.019 


9.263 


0.023 


2MASS H 


10.124 


0.023 


8.860 


0.025 


2MASS K, 


9.835 


0.021 


8.598 


0.025 


GL97 J 


10.68 


0.05 






GL97H 


10.12 


0.05 






GL97K 


9.79 


0.04 






VdSOO J 






9.27 


0.05 


VdSOO H 






8.81 


0.05 


VdSOO K 






8.61 


0.05 


VdSOO L 






8.27 


0.10 




F„ 


AF V 


F„ 


AF„ 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


MSX A 8.28 






0.72 


0.03 


MSXC 12.13 






3.35 


0.18 


MSX D 14.65 






7.1 


0.4 


MSX E 21.34 






17.7 


1.1 


IRAS 12 


2.23 


0.11 


3.31 


0.20 


IRAS 25 


9.83 


0.59 


30.6 


1.2 


IRAS 60 


3.62 


0.36 


17.1 


1.9 



Table 6. The total extinction A v and the bolometric luminosity L of the 
program stars. 



IRAS 


A v 


L at D = 5 kpc 




mag 


L G 


08281-4850 


4.33 


2450 


14325-6428 


3.32 


7250 




-0.5 0.0 0.5 1.0 1.5 

log(X) (urn) 



Fig. 2. The spectral energy distribution of IRAS 14325-6428. 
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Fig. 1. The spectral energy distribution of IRAS 08281-4850. 



Fig. 3. The spectra of IRAS 14325-6428 and IRAS 08281-4850 com- 
pared with the A7Iab standard HD 81471 and the non-enriched post- 
AGB star HD 133656. 



Reyniers & Van Winckel 2003; Van de Steene & van Hoof 
2003). IRAS 14325-6428 is observed with UVES on the VLT 
UT2 telescope (Kueyen), as a member of a larger sample of 
seven post-AGB objects that were observed in service mode dur- 
ing ESO period #73. IRAS 08281-4850 was observed in visitor 
mode with the EMMI echelle grating #9 with cross disperser 
#3 at the NTT during period #70. The resolving power of the 
UVES spectra varies between ~55,000 and ~60,000. The resolv- 
ing power of our EMMI spectra is significantly lower (~8,400). 
The spectral interval covered and some other details about the 
observations are given in Table 7. 
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Table 7. Log of the high-resolution observations. For IRAS 08281- 
4850, a spectral gap occurs between 493 nm and 500 nm. For 
IRAS 14325-6428, spectral gaps occur between 577 nm and 583 nm and 
between 854.4 nm and 864.5 nm due to the spatial gap between the two 
UVESCCDs. 



date UT exp.time wavelength S/N 

start (sec) interval (nm) 
IRAS 08281-4850 (NTT + EMMI) 

2003- 02-04 01:11 9000 395-795 75 
IRAS 14325-6428 (VLT-UT2 + UVES) 

2004- 05-13 07:13 1800 374.5-498 130 
2004-05-13 07:13 1800 670.5-1055 170 
2004-05-13 06:39 1800 477.5-681 190 



The reduction of our UVES spectra was performed in the 
dedicated "UVES context" of the midas environment and in- 
cluded bias correction, cosmic hit correction, flat-fielding, back- 
ground correction and sky correction. We used average extrac- 
tion to convert frames from pixel-pixel to pixel-order space. The 
reduction of our EMMI spectrum was done within the echelle 
package in iraf following the user's guide by D. Willmarth and 
J. Barnes (1994). The spectra were normalised by dividing the 
individual orders by a smoothed spline function defined through 
interactively identified continuum points. For a detailed descrip- 
tion of the reduction procedure, we refer to Reyniers (2002). In 
Table 7, we also list some indicative signal-to-noise values of the 
final data product. 

Sample spectra of our programme stars can be found in 
Figs. 3 and 6. In Fig. 3, the spectra of IRAS 08281-4850 
and IRAS 14325-6428 are compared with the spectra of the 
A7Iab standard HD 81471 and the non-enriched post-AGB 
star HD 133656. The spectrum of HD 81471 is retrieved from 
UVESPOR the library of high-resolution UVES spectra of stars 
across the Hertzsprung-Russell diagram (Bagnulo et al. 2003). 
HD 133656 is discussed in Van Winckel et al. (1996). The spec- 
trum of HD 133656 is an ESOL5-m+FEROS spectrum taken on 
June 27, 2001. This object has atmospheric parameters com- 
parable with the program stars (T e ff, \ogg, = (8000 K, 1.0, 
3.0 kms -1 ) but is slightly more metal deficient ([Fe/H] = -0.9). 
The s-process enhancement of the two IRAS stars is clear from 
the Ban line. IRAS 08281-4850 is the strongest s-enriched one, 
but the stronger lines are also due to a slightly lower temperature 
of this object. 

4. Abundance analysis and results 

4.1. Method 

The general methodology is already extensively discussed in our 
previous papers (e.g. Reyniers et al. 2004; Deroo et al. 2005) and 
will not be repeated here. Here we only remind that we make use 
of the latest Atlas models (Castelli & Kurucz 2004) in combina- 
tion with the latest version (April 2002) of Sneden's LTE line 
analysis program MOOG (Sneden 1973). The value of the mi- 
croturbulent velocity of IRAS 14325-6428 derived from the Fe n 
lines is high (£ t = 10km s~'), but this high value is also supported 
by the results of other ions as can be seen in Fig. 4. 

Since the spectral resolution of the UVES spectra is much 
higher than that of the EMMI spectrum, we started the analy- 
sis with IRAS 14325-6428. In order to avoid undetected blends 
that are caused by the lower spectral resolution of the EMMI 
spectrum, we took the line list of IRAS 14325-6428 as our start - 
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Fig. 4. The relatively high microturbulent velocity (£ = 10 km s -1 ) for 
IRAS 14325-6428 is also supported by the results of ions other than 
Fen. 



ing point for the analysis of IRAS 08281-4850. As a conse- 
quence, all lines that we used for IRAS 08281-4850 are also 
present in the line list of IRAS 14325-6428. This choice en- 
sures the consistency of the analyses, implying that a compar- 
ison between the abundance results is highly reliable. Due to 
the stronger s-process enrichment of IRAS 08281-4850 com- 
pared to IRAS 14325-6428, the spectrum of IRAS 08281-4850 
shows significantly more lines of s-process elements. Therefore, 
we performed an extra search for s-process lines in the spectrum 
of IRAS 08281-4850. This explains why we have used 12 Lan 
lines and 1 Ndn in the analysis of IRAS 08281-4850, against 
only 5 La n lines and no Nd n lines for IRAS 14325-6428. 

4.2. Abundance results 

The final results of our stellar atmosphere parameter determina- 
tion and our abundance analysis are listed in Table 8. The first 
column of this table gives the actual ion. Then, for each star, 
the following columns are listed: the number of lines used; the 
mean equivalent width in mA; the absolute abundances by num- 
ber loge = log(N(el)/N(H))+12; the line-to-line scatter cr\ t \; the 
abundance relative to iron [el/Fe], and an estimate of the total 
uncertainty on the abundance <x tot (see Sect. 4.3). For the ref- 
erences of the solar abundances (the middle column) needed to 
calculate the [el/Fe] values: see Reyniers et al. (2007). 

The abundances are also graphically presented in Fig. 5. On 
this figure, the different groups of elements are marked with dif- 
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Table 8. Abundance results for IRAS 08281-4850 and IRAS 14325-6428. For the explanation of the columns: see Sect. 4.2. 
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IRAS 08281-4850 

T eff = 7750K 
log g= 1.0 (cgs) 
£ = 4.5 kms" 1 



IRAS 14325-6428 

T eff = 8000K 
logg= 1.0 (cgs) 
£= 10.0 km s 1 



ion 


N 


W A 


loge 


cm 


[el/Fe] 


0"tot 


sun 


N 


W A 


log 6 


0" M 


[el/Fe] 


O-fot 


Ci 


12 


117 


9.17 


0.24 


0.93 


0.15 


8.57 


30 


64 


9.09 


0.12 


1.07 


0.13 


Ni 


2 


87 


8.21 


0.14 


0.55 


0.20 


7.99 


6 


94 


7.98 


0.18 


0.54 


0.14 


Oi 


2 


84 


8.89 


0.10 


0.36 


0.22 


8.86 


10 


72 


8.89 


0.11 


0.58 


0.14 


Nai 


1 


61 


6.76 




0.76 


0.25 


6.33 


2 


16 


6.50 


0.02 


0.72 


0.20 


Mgi 


1 


107 


7.31 




0.10 


0.24 


7.54 


3 


72 


7.52 


0.13 


0.53 


0.18 


Mgn 














7.54 


3 


56 


7.25 


0.09 


0.26 


0.19 


Ah 














6.47 


1 


147 


5.74 




-0.18 


0.23 


Si 














7.33 


1 


17 


7.26 




0.48 


0.24 


Cai 


2 


72 


6.40 


0.08 


0.37 


0.20 


6.36 


2 


18 


6.20 


0.02 


0.39 


0.20 


Sen 


3 


110 


3.25 


0.11 


0.41 


0.12 


3.17 


4 


55 


3.01 


0.07 


0.39 


0.11 


Tin 


9 


113 


4.81 


0.15 


0.12 


0.07 


5.02 


26 


84 


4.58 


0.11 


0.11 


0.04 


Cri 


1 


122 


5.60 




0.26 


0.22 


5.67 


2 


61 


5.34 


0.10 


0.22 


0.16 


Cm 


9 


70 


5.40 


0.26 


0.06 


0.10 


5.67 


17 


58 


5.23 


0.11 


0.11 


0.04 


Mnn 














5.39 


2 


24 


5.07 


0.11 


0.23 


0.15 


Fei 


19 


79 


7.29 


0.14 


0.11 


0.11 


7.51 


34 


37 


7.07 


0.13 


0.11 


0.11 


Fen 


15 


98 


7.18 


0.11 


0.00 




7.51 


22 


74 


6.96 


0.11 


0.00 




Nin 














6.25 


3 


33 


5.67 


0.24 


-0.03 


0.12 


Yn 


9 


116 


3.81 


0.17 


1.90 


0.08 


2.24 


15 


70 


3.01 


0.19 


1.32 


0.06 


Zrn 


1 


170 


3.85 




1.58 


0.26 


2.60 


14 


91 


3.19 


0.14 


1.14 


0.06 


Ban 














2.13 


1 


84 


3.13 




1.55 


0.20 


La ii 


12 


75 


2.93 


0.18 


2.13 


0.08 


1.13 


5 


54 


1.89 


0.10 


1.31 


0.06 


Cen 














1.58 


2 


16 


2.15 


0.05 


1.12 


0.15 


Ndn 


1 


52 


2.78 




1.61 


0.21 


1.50 














Ndm 


1 


101 


2.57 




1.40 


0.22 


1.50 


3 


60 


2.00 


0.11 


1.05 


0.14 


Smn 














1.01 


1 


11 


1.69 




1.23 


0.20 



summary 


summary 


[Fe/H] = -0.33 


[Fe/H] = -0.55 


C/0=1.9 


C/0= 1.6 


[a/Fs] = +0.2 (a: Mg, Ca, Ti) 


[a/Fe] = +0.3 (a: Mg, S, Ca, Ti) 


[ls/Fe] = +1.74 


[ls/Fe] = +1.23 


[hs/Fe] = +1.91 


[hs/Fe] = +1.29 


[hs/ls] = +0.17 


[hs/ls] = +0.06 



ferent symbols. We will summarize the main results for each of 
these groups. 

Metallicity Both stars are mildly metal deficient, with iron 
abundances of [Fe/H] = -0.33 and -0.55 for IRAS 08281- 
4850 and IRAS 14325-6428 respectively. The other iron 
peak elements, if present, follow this deficiency. 

CNO-elements Both stars are clearly carbon enriched, with an 
enrichment around [C/Fe] +1 for both stars. As a conse- 
quence, we derive also high C/O number ratios for both stars. 
One has to note, however, that in IRAS 08281-4850 the un- 
certainty on the oxygen abundance prevents an accurate C/O 
number ratio for this star. 

a-elements The simple mean of the [el/Fe] values of the 
(available) a-elements yields [a/Fe]=+0.2 and +0.3 for 
IRAS 08281-4850 and IRAS 14325-6428 respectively. Such 
an enhancement is normal for stars in this metallicity range, 
as a consequence of the galactic chemical evolution and 
therefore does not correspond to an intrinsic enhancement. 

s-process elements It is clear that the s-process enrichment of 
the two objects under study is very strong. 

The s-process elements observed in evolved stars can be divided 
into two groups: the light s-process elements around the magic 
neutron number 50 (Sr, Y, Zr) and the heavy s-process elements 



around the magic neutron number 82 (Ba, La, Ce, Pr, Nd, Sm). 
Three s-process indices are generally defined: [ls/Fe], [hs/Fe] 
and [hs/ls]. To be consistent with our earlier papers on simi- 
lar stars (Van Winckel & Reyniers 2000; Reyniers et al. 2004; 
Reyniers & Cuypers 2005), we define the ls-index as the mean 
of the Y and Zr abundances and the hs-index as the mean of the 
Ba, La, Nd and Sm abundances, with unavailable elements esti- 
mated using the tables of Malaney (1987). All indices are listed 
in Table 8. 



4.3. Uncertainty estimates 

For the error analysis, we followed the same method as described 
in Deroo et al. (2005). We slightly changed formula (1) of this 
paper, in the sense that for the uncertainty induced by the model 
Cmod, we confined the parameter space to consistent models, i.e. 
models for which there is ionisation equilibrium between Fe i 
and Fe n. To be more precise, in order to calculate <x mo( j, we stud- 
ied the abundance changes for two different consistent models 
(T eff =7750K, logg=0.6) and (T eff =8250K, log £=1.5), together 
with a change of the microturbulent velocity of £,=2 km s _1 . The 
total uncertainty on the [el/Fe] abundances <x tot can be found in 
Table 8, and is the quadratic sum of the uncertainty on the mean 
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Fig. 5. [el/Fe] values for IRAS 14325-6428 (upper panel) and Fig. 6. Four clear diffuse interstellar bands in the spectrum of 
IRAS 08281-4850 (lower panel). The uncertainty on the [el/Fe] values IRAS 14325-6428. A spectral synthesis is overplotted in red, to facil- 
is the total uncertainty cr tot as listed in Table 8. itate the DIB-detection. 



due to line-to-line scatter, the uncertainty induced by the model, 
and the uncertainty on the Fe abundance: 

/. 0"ltl " ~~ cr Fe 

cr tot = /(—= ) 2 + (o- mod ) 2 + (^==) 2 

If less than 5 lines were available, a line-to- line scatter of 0.2 dex 
was applied. 

5. Diffuse Interstellar Bands 

Diffuse Interstellar Bands (DIBs) are broad absorption lines of 
interstellar origin that are seen in the spectra of reddened objects. 
The carriers of these DIBs are still not known, but polycyclic 
aromatic hydrocarbons (PAHs) are amongst the most probable 
candidate carriers. Due to the severe mass loss in the preceding 
AGB phase, post-AGB stars are often enshrouded by carbon-rich 
circumstellar dust, causing severe reddening. Therefore, post- 
AGB stars are ideal testlabs to search for possible circumstel- 
lar DIBs. If, on the other hand, the DIBs are detected to be in- 
terstellar, a rough division can be made between the interstellar 
and circumstellar component of the total reddening towards the 
post-AGB star, since some DIBs correlate quite well with the 
interstellar reddening. 

5.1. DIBs in IRAS 14325-6428 

During our analysis of IRAS 14325-6428, we realised that strong 
DIBs are indeed present in this spectrum. We initiated a system- 
atic search for DIBs. First, we selected eight well-known and 
widely studied DIBs from the list in Herbig (1995). For each of 
these DIBs, we made a spectral synthesis in the vicinity of the 
DIB wavelength, based on the abundances found in our abun- 
dance analysis (Table 8). The DIB can then easily be identified as 
the spectral feature that is not fitted. Four DIBs are shown in Fig. 

6. The studied DIBs in the spectrum of IRAS 14325-6428 are 



listed in Table 9. The columns of this table represent: rest wave- 
length of the DIB; the observed wavelength in IRAS 14325- 
6428; the heliocentric radial velocity of the DIB; the measured 
equivalent width; the intrinsic strength of the DIB as seen in the 
DIB standard HD 183143; and the (interstellar) reddening de- 
rived from the latter relation. The reference for columns (1) and 
(5) is Herbig (1995). 

The heliocentric radial velocities of the DIBs show some 
scatter, but they are definitely different from the heliocentric 
velocity of the star (-87 kms _I ) and show velocities very dif- 
ferent from any expected outflow of the circumstellar mate- 
rial. Therefore, the DIBs seen in the spectrum of IRAS 14325- 
6428 are the imprint of an interstellar cloud between the object 
and the observer. A rough estimate for the distance towards the 
DIB producing cloud is obtained through the formula of Lang 
(1980), which describes the rotation of the Galactic Plane, yield- 
ing 1.6 kpc. Obviously, this distance is also a lower limit for the 
distance towards IRAS 14325-6428. The scatter in the DIB ve- 
locities could be caused by different DIB producing clouds in the 
line of sight towards IRAS 14325-6428, but also the rest wave- 
lengths of the DIBs are difficult to quantify due to their complex 
fine structure (e.g. Galazutdinov et al. 2003). 

For most DIBs, there is a correlation between its strength 
and the reddening that is caused by the DIB producing cloud. 
The DIB standard that is often used to quantify this relation, is 
HD 183143 (E(B-V)= 1.28). In Table 9, the DIB strengths of 
this standard, normalized to E(B-V)= 1, are given in column 
5, together with the inferred reddening for the DIB producing 
cloud towards IRAS 14325-6428 (column 6). A mean interstellar 
reddening of E(B-V) = 0.6 is obtained. If we combine this with 
the total reddening towards IRAS 14325-6428 of E(B-V) tot =1.1 
(Table 6), we can conclude that a significant fraction (around 
50 %) of the total reddening towards IRAS 14325-6428 is of in- 
terstellar origin. 
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Table 9. Studied DIBs in the spectrum of IRAS 14325-6428. 



DIB A 


obs. A 


Vhelio 


EW 


EW 

E(B-V) 


E(B-V) 


(A) 


(A) 


(kms- 1 ) 


(A) 


HD 183143 


from DIB 


5780.45 


5779.43 


-16.6 


0.417 


0.626 


0.7 


5796.98 


5796.15 


-6.3 


0.087 


0.186 


0.5 


5849.65 


5849.57 


-1.1 


0.042 


0.064 


0.7 


6195.95 


6195.61 


-13.2 


0.043 


0.063 


0.7 


6283.86 


6283.60 


-9.2 


1.043 


1.520 


0.7 


6379.20 


6378.95 


-8.5 


0.041 


0.096 


0.4 


6613.62 


6613.28 


-12.4 


0.138 


0.280 


0.5 


6993.07 


6992.82 


-7.7 


0.079 


0.142 


0.6 


mean (±<x) 




-9.4 (±4.8) 






0.6 (±0.1) 



5.2. DIBs in IRAS 08281 -4850 

Due to the lower resolution and the lower S/N ratio of the 
IRAS 0828 1 -4850 spectra, we were not able to do a similar study 
for this source. A qualitative comparison of the strength of the 
strong DIB at 6 196 A indicates that IRAS 08281-4850 has cer- 
tainly weaker DIBs than IRAS 14325-6428. The total reddening 
towards IRAS 08281-4850 is, however, stronger than the one to- 
wards IRAS 14325-6428, implying that the circumstellar com- 
ponent of the total reddening of IRAS 08281-4850 is probably 
large. 



6. Discussion 

IRAS 08281-4850 and IRAS 14325-6428 were selected on the 
basis of their position in the IRAS colour-colour diagram 
(Pottasch et al. 1988) and the lack of free-free radio continuum 
emission (Van de Steene & Pottasch 1993, 1995). Both objects 
remained poorly studied until now. Our chemical analysis pre- 
sented in this paper shows that both objects are among the hottest 
members of the s-process enhanced post-AGB stars known to 
date and illustrates that the rich IRAS legacy of objects in the 
transition between the AGB and the PN phase is far from being 
harvested. Although the distance to the objects is not well con- 
strained, both objects do not show an extremely high luminosity 
(Table 6) and, given the sub-solar metallicity, both objects repre- 
sent the final evolutionary phase of a star with a low initial mass 
(Mi<2M ). 

There is general agreement that the source of the neutrons 
for the s-process in low and intermediate AGB stars is the " 13 C 
source". A long standing problem is the formation of the 13 C 
pocket itself. Stellar models using a standard treatment of mixing 
cannot reproduce the 13 C pocket at a level which is high enough 
for the s-process to take place (Herwig 2005). Until now, there is 
no satisfying description for this phenomenon, and the different 
modeling groups use different descriptions. Many models use 
some kind of overshoot mechanism (e.g. Herwig et al. 1997). 
Also differential rotation has recently been studied as a possible 
mixing mechanism driving the proton engulfment, but this rota- 
tionally induced mixing alone cannot account for the formation 
of a large enough 13 C pocket (Herwig et al. 2003). Other groups 
avoid this formation problem by assuming an ad-hoc 13 C pocket 
in the He intershell (models by e.g. Gallino et al. 1998) or an ad- 
hoc proton density profile (models by e.g. Goriely & Mowlavi 
2000). The observed spread in s-process efficiency (see Sect. 1) 
is then reproduced by a variable 13 C pocket strength. Since the 
13 C pocket is build upon primary synthesised 12 C, the formation 
is thought to be largely independent of the initial metallicity and 
a richer s-process nucleosynthesis is expected for lower metallic- 



ity stars, because more neutrons become available per iron-seed. 
In practise, such a trend has not been found for post-AGB stars. 

The very recent results of the population synthesis mod- 
els by Bonacic Marinovic et al. (2006) and Bonacic Marinovic 
et al. (2007) are very interesting in this context. Their approach 
is to combine stellar population synthesis with a rapid stellar 
evolution code including AGB nucleosynthesis and evolution. 
Interestingly, their models do reproduce both the observed di- 
chotomy (see Fig. Id in Bonacic Marinovic et al. 2006), and the 
observed spread. Contrary to previous models, only a limited 
spread in the strength of the 13 C pocket is needed to reproduce 
the observed efficiency spread. With the new results presented 
in this paper, we are able to add two important datapoints to 
constrain this promising new generation of models. Indeed, the 
two stars discussed here show a strong s-process enrichment, 
although they are only mildly metal deficient. Particularly the 
results of IRAS 08281-4850 are on the predicted boundaries of 
both the [hs/ls] index (Fig. 3 of Bonacic Marinovic et al. 2007) 
and the Zr overabundance (Fig. 8 of Bonacic Marinovic et al. 
2007). 

7. Conclusion 

With the comprehensive abundance analysis presented in this 
paper, the objects IRAS 08281-4850 and IRAS 14325-6428 join 
the group of post-AGB stars with a clear post third dredge-up 
signature, since they do not only show a clear carbon enhance- 
ment, but also a strong enrichment in s-process elements. This 
enrichment is surprisingly strong with a high ratio of heavy ver- 
sus light s-process elements, despite the only mild metal defi- 
ciency of both objects. This provides additional evidence of the 
intrinsic s-process efficiency spread in (post-)AGB stars. A sys- 
tematic analysis of post-AGB stars in the PNe locus of the IRAS 
colour-colour diagram would be a very rewarding program to 
study the relation between the AGB nucleosynthesis, the dredge- 
up efficiency and the overall stellar evolution of the central star. 
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Table 4. Absolute flux calibration for the various photometric bands 
used in this paper. 
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Table 5. The spectral energy distribution (sorted by wavelength) given 
as AF A in SI units. 
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